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SEQOPT 

A  SOLUTION  TO  THE  FORCE  STRUCTURING  PROBLEM 

ABSTRACT 

SEQOPT  has  been  designed  to  solve  a  specific  force  structuring 
and  allocation  problem.  Given  the  following: 

1.  A  set  of  targets  each  with  a  relative  value  and  an  asso¬ 
ciated  single  shot  kill  probability; 

2.  A  set  of  launch  points  for  weapons  and  the  number  of 
weapons  allowed  at  each  launch  point; 

3.  The  payload  types  available  for  use  in  these  weapons; 

4.  The  set  of  accessible  targets  for  a  weapon  carrying  each 
payload  type  from  every  launch  point. 

SEQOPT  will  then  determine  the  payload  type  for  each  weapon 
in  the  force,  and  the  allocation  of  each  weapon  to  its  set  of  accessible 
targets,  in  a  way  which  attempts  to  maximize  the  expected  damage  to  the 
target  complex. 
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SEQOPT 

A  SOLUTION  TO  THE  FORCE  STRUCTURING  PROBLEM 

1.  INTRODUCTION 

SEQOPT  is  a  computer  model  written  in  FORTRAN  IV  for  the 
BRLESC  II  (Ballistic  Research  Laboratory  Electronic  Scientific  Computer) 
which  chooses  the  payload  type  to  be  used  by  each  weapon  in  a  force,  and 
allocates  these  weapons  to  a  target  complex  in  a  manner  whose  aim  is  to 
maximize  the  expected  damage  of  a  target  complex. 

The  method  used  in  the  model  is  called  sequential  optimization 
and  does  not  necessarily  yield  a  globally  optimal  solution.  Those 
methods  which  would  yield  an  optimal  solution  are  at  this  point  in  time 
unfeasible  for  computer  use. 


2.  DEFINITIONS 


Before  the  problem  can  be  adequately  described,  some  terms, 
which  will  have  precise  meanings  in  this  report,  must  be  defined: 

Relative  Target  Value:  Each  target  has  associated  with  it  a 
relative  value,  which  can  be  based  upon  any  consistent  attri¬ 
bute,  such  as  industrial  floor  space,  population,  etc.  (to  be 
determined  by  the  analyst) . 

Weapon :  The  weapon  can  be  assumed  to  be  a  bomber  or  a  varia¬ 
ble  payload  missile  which  is  capable  of  depositing  the  war¬ 
heads  it  carries  on  one  or  more  of  the  targets  which  are 
accessible  to  it. 


Payload  Type:  The  payload  is  composed  of  a  specific  number 
of  warheads,  all  of  which  are  exactly  alike  in  characteristics 
Thus  payload  types  differ  not  in  the  types  of  warheads  in¬ 
volved  but  in  the  number  of  warheads  involved. 


Warhead:  The  warhead  could  be  either  a  bomb  or  an  atomic 
warhead,  specified  by  yield;  the  standard  warhead  used  in  all 
payload  types . 


Single  Shot  Kill  Probability  (PK) : 

each  target  is  actually  a  function 
warhead  characteristics.  However, 
fixed  throughout  the  problem,  a  PK 
calculated. 


The  PK  associated  with 
of  target  and  attacking 
since  the  warhead  type  is 
for  each  target  can  be 
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Launch  Point:  The  launch  point  can  be  considered  to  be  a 
bomber  base  or  the  point  from  which  missiles  are  launched 
from  a  submarine.  Capacity  refers  to  the  number  of  weapons 
which  are  to  be  launched  from  a  launch  point. 

Accessible  Targets:  Associated  with  each  weapon  carrying  a 
specific  payload  type  is  a  maximum  range.  The  accessible  tar¬ 
gets  for  a  weapon  with  a  specific  payload  from  a  given  launch 
point  are  all  those  targets  in  the  target  complex  lying  within 
the  maximum  and  minimum  ranges  associated  with  that  payload 
type. 

Out load:  The  outload  is  the  choice  of  payload  type  for  each 
weapon  at  each  launch  point,  the  final  structure  of  the  force. 


3.  MISSILE  ALLOCATION  PROBLEM 
3. 1  Description  of  Problem. 

A  great  deal  of  work  has  been  done  in  solving  the  so  called 
"Missile-Allocation  Problem."  The  problem  considered  is  this:  given 
an  existing  weapon  force  and  set  of  targets,  what  is  the  optimal  alloca¬ 
tion  of  weapons  to  targets?^  (i.e..  How  can  we  maximize  expected  damage 
to  the  target  complex?) 

The  SEQOPT  model  is  concerned  with  a  variation  of  this  stand¬ 
ard  allocation  problem.  The  new  problem  may  be  described  as  follows: 
given  a  set  of  fixed  launch  points,  the  number  of  weapons  at  each  launch 
point,  and  a  target  complex,  choose  the  payload  type  for  each  weapon, 
and  allocate  the  weapons  to  targets  so  as  to  maximize  the  expected 
damage  to  the  target  complex. 

The  obvious  difference  between  these  problems  is  that,  in  the 
standard  allocation  problem,  we  are  given  a  fixed  force  and  asked  to 
allocate  it  optimally.  In  the  force  structuring  problem,  which  SEQOPT 
solves,  we  are  asked  to  actually  structure  the  force,  so  that  when  it  is 
allocated  the  expected  damage  is  optimal. 


^Matlin,  Samuel,  "A  Review  of  the  Literature  on  the  Missile  Allocation 
Problem,"  Operations  Research  18,  pp.  334-373  (1970). 
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An  example  will  serve  to  clarify  the  distinction  between  prob¬ 
lems,  and  what  is  meant  by  "structuring  the  force."  Suppose  the  offense 
consists  of  one  launch  point  which  has  a  capacity  for  two  weapons 
(bombers  or  missiles) : 

1.  The  standard  problem  would  state  that  these  weapons  carry 
payloads  consisting  of  N^,  and  N2  warheads,  respectively.  Given  the  tar¬ 
gets  which  are  accessible  to  weapons  carrying  these  payloads,  allocate 
these  weapons  to  the  target  complex  optimally. 

2.  The  new  problem  asks:  If  the  payloads  for  the  two  weapons 
can  vary  between  A  and  B  warheads  per  weapon,  determine  the  payloads  for 
the  two  weapons  which,  when  allocated  optimally,  will  yield  the  greatest 
expected  damage  to  the  target  complex.  Notice  here  that  and  N2,  as 
well  as  their  allocations,  are  determined.  (A  <_  N^,  N2  —  ^  * 

Obviously,  the  expected  value  damaged  achieved  in  the  new 
problem  is  at  least  as  great  as  that  in  the  standard  problem.  This  is 
because,  in  the  new  problem,  all  possible  combinations  of  payloads  are 
tested  for  return;  then  the  outload  which  yields  the  greatest  damage  is 
chosen. 

In  general,  it  will  be  found  that  the  maximum  return  is  not 
achieved  by  carrying  the  "heaviest"  payloads  possible.  This  is  due  to 
the  fact  that,  as  the  payload  size  increases,  the  maximum  effective 
range  of  the  weapon  decreases  under  the  condition  of  a  constant  quantity 
of  propellant.  Therefore,  when  carrying  the  maximum  payload,  some  of 
the  most  valuable  targets  may  lie  out  of  range;  consequently,  some 
value  is  not  damaged. 

The  exception  to  this  rule  exists  when  all  targets  are  within 
range  of  weapons  carrying  the  heaviest  payloads  from  every  launch  point 
considered.  In  this  case  it  is  advantageous  to  use  only  the  heaviest 
payloads  available,  since  the  use  of  lighter  payloads  makes  no  new  tar¬ 
gets  available.  However,  this  special  case  is  of  little  interest  and 
so  we  will  not  consider  it. 


9 


Precise  mathematical  descriptions  of  the  standard  allocation 
and  force  structuring  problems  follow. 

3.2  Mathematical  Description  of  Allocation  Problems. 

Some  definitions  are  needed  before  we  can  formulate  the 

problem: 

LAUNCH  POINT  VARIABLES 

W.  =  the  number  of  weapons  allowed  at  the  i*'*1  launch  point 
(capacity)  . 

LP  =  the  number  of  launch  points  being  considered. 

6.  .  =  the  number  of  weapons  of  the  payload  type  placed  at  the 

. th  .  ,  . 

i  launch  point. 

A.  .,  =  1  if  the  target  is  in  range  of  the  j™  payload  type 
1JK  th 

from  the  i  launch  point;  0  otherwise. 

A 

A^  =  the  vector  whose  components  are  A„^  f°r  k=l,...,T, 


Aij  =  {Aijl»  Aij  2  *  *  *  * ,AijT}’ 


WEAPON  VARIABLES 

M  =  the  number  of  payload  types  being  considered, 

cu  =  the  number  of  warheads  in  the  j  payload  type. 

TARGET  VARIABLES 


X.  .. 
ljk 


the  number  of  warheads  from  weapons  of  the  j*"*1  payload  type 

th  th 

from  the  i  launch  point,  placed  on  the  k  target. 

the  number  of  warheads  placed  on  the  k^  target 


LP 


M 

Z 


N.  =  Z  Z  X.  .. 
k  ......  ijk 

1=1  j=l 


=  the  number  of  targets  in  target  complex. 


th 


f,  =  the  expected  damage  function  associated  with  the  k  target 

PK^  =  the  single  shot  kill  probability  of  the  k  target  (based 
on  target  and  weapon  characteristics) . 

V,  =  the  relative  value  of  the  k^1  target. 
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Both  allocation  problems  have  the  same  formulation,  the 
same  objective  function,  but  a  slight  variation  in  the  constraints. 

T 

Maximize  G  =  I  f^V.,  P!^,  H  ) 
k=l 

Subject  to: 


111 

+  X112  + 

. . .  +  x11T  <$n  , 

, 

121 

+  X122  + 

. . .  +  x12T  =  <$12  , 

,  a2 

+  X 


klM2 


+  X 


1MT 


=  <5 


lm  '  °M 


where  6, .  =  0,1,, 
lj 

M 

and  Z  <S.  .  =  W, 
j=i  »  1 


LP11  + 

X  + 

LP12  *• 

.  +  XLpiT  6Lpi  . 

°1 

LP21  + 

X  + 

LP22 

’  +  XLP2T  =  <SLP2  * 

“2 

where  <SLp_.  =  0,1,2,. .. ,WLp 
M 

LPM1  + 

X  + 

LPM2 

’  +  XLPMT  =  6LPM  * 

“M 

md  ,E,  SLPj  ’  WLP 

J-l 

All  X...  >  0,  and  if  A...  =  0,  then  X...  =  0. 
ljk  —  ijk  ’  ljk 


The  entire  difference  between  the  two  problems  lies  in  the  6^  variables, 

which  determine  the  outload  or  structure  of  the  force.  In  the  old  problem 

the  are  all  fixed  at  the  outset  of  the  problem,  and  they  satisfy  all 

constraints  associated  with  them.  The  new  problem  asks  the  analyst  to 

determine  the  best  choice  of  6..,  i.e.,  the  force  which  when  allocated 

il 

to  the  targets  optimally,  will  maximize  expected  damage. 

Thus,  the  allocation  problem  asks  only  for  the  set  of  |x.^j 
which  maximize  expected  value  damaged  given  ^6^  j.  The  force  structur¬ 
ing  problem  asks  the  analyst  to  determine  and  ^Xijk)  max*- 

mize  the  same  objective  function. 
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It  may  be  that  there  are  not  enough  warheads  available  to  fill 
the  force  determined.  The  analyst  then  has  two  choices;  either  lower  the 
number  of  weapons  to  be  used,  or  lower  the  number  of  warheads  in  the  pay- 
load  types. 

3. 3  Discussion  of  Constraints. 

First  notice  that  there  is  a  set  of  equations  present  for  each 
launch  point.  These  ensure  that  the  number  of  warheads  from  weapons  of 
a  given  type,  allocated  to  its  accessible  targets,  does  not  exceed  the 
number  of  warheads  in  weapons  of  that  type  present  at  that  launch  point 


The  6^.  are  restricted  to  integers  since  they  represent 
weapons  of  a  specific  type.  The  restriction  that 


M 

Z  6 . .  =  W.  for  all  i 

j-1  1J  1 


ensures  that  the  capacity  of  launch  point  i  is  not  exceeded. 


Finally,  the  fact  that  if  =  0,  then  roust  also  equal 

0,  is  obvious  when  one  considers  that  a  weapon  cannot  attack  a  target 
lying  outside  its  maximum  range. 


There  is  one  real-world  constraint  which  is  not  considered 

in  this  model.  This  is  concerned  with  the  set  of  targets  which  are 

within  range  of  a  certain  payload  type  from  a  fixed  launch  point;  that 

is  all  K.  such  that  A...  =1  for  fixed  i.j.  The  model  assumes  that  a 
3  ljk 


th  th 

weapon  carrying  the  j  payload  type  from  the  iu  launch  point  can 
attack  any  subset  of  its  accessible  targets  (subject  of  course  to  the 
constraint  that  it  cannot  attack  more  targets  than  warheads  it  carries) . 
This  simplifying  assumption  is  in  general  not  true,  because  the  number 
of  targets  that  can  be  attacked  depends  upon  such  things  as  time  of 
flight  to  the  first  and  subsequent  targets  and  whether  fuel  restrictions 
are  exceeded.  The  problem  concerned  with  the  sequencing  of  targets  to 
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be  attacked  so  as  to  remain  within  fuel  and  time  of  flight  constraints 
is  known  as  footprinting. 

Once  a  footprint,  or  a  set  of  targets  to  be  attacked  by  one 

2 

weapon  has  been  determined,  the  Branch  and  Bound  Algorithm,  developed 
by  Little,  et  al.,  can  be  used  to  determine  the  least  distance  route 
through  all  targets  in  the  set.  The  results  of  this  algorithm  when  com¬ 
pared  with  maximum  time  of  flight  can  be  used  to  determine  the  achieva- 
bility  of  a  given  footprint. 

3 . 4  Objective  Function. 

Note  that  the  objective  function  G  is  a  sum  consisting  of  one 
expected  damage  function,  f^.,  associated  with  each  target.  Each  f^,  is 
a  function  of  three  variables: 

1.  The  initial  relative  value  of  the  k^1  target,  V,  . 

th  ^ 

2.  The  single  shot  kill  probability  of  the  k  target,  PK,  . 

th  ^ 

3.  The  number  of  warheads  applied  to  the  k  target,  N^,. 

There  are  various  forms  of  these  expected  damage  functions; 
however,  the  following  functional  form  is  used  in  SEQOPT : 

A  ■  l1  -  (1-pKk)NK]  •  \ 

These  functions  are  all  monatomically  increasing,  bounded 
above  by  V^.,  and  follow  the  law  of  diminishing  returns. 


EXPECTED 

DAMAGE 


NUMBER  OF  WARHEADS  APPLIED 

2Little,  J.  D.  C.,  et  al . ,  "An  Algorithm  for  the  Traveling  Salesman 
Problem,"  Operations  Research  11,  pp.  972-989  (1963). 
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3.5  Possible  Methods  of  Solution  to  Allocation  Problems. 

The  body  of  mathematical  programming  techniques  was  first 
explored  for  a  means  of  solving  the  force  structuring  problem.  Linear 
programming  was  immediately  eliminated  due  to  the  non-linearity  of  the 
objective  function. 

Next,  dynamic  programming  was  considered,  and  a  solution  to 
the  force  structure  problem  was  found.  Embedded  within  the  recursive 
relationship  of  this  solution  is  the  dynamic  programming  solution  of  the 
standard  allocation  problem. 

Thus,  in  theory,  dynamic  programming  can  be  used  to  solve  the 
force  structuring  problem;  however,  the  memory  requirements  are  so 
large  that  the  method  is  at  present  unfeasible.  (The  interested  reader 
should  consult  the  appendix  of  this  report  for  the  dynamic  programming 
formulation  which  solves  these  problems.) 

Next  the  Branch  and  Bound  Algorithm,  developed  by  Little, 

2 

et  al.,  to  solve  the  "Traveling  Salesman  Problem"  was  considered. 

This  method  divides  the  possible  set  of  into  a  tree- like  structure. 
If  a  bound  could  be  placed  on  these  branches,  large  numbers  of  possible 
solutions  could  be  eliminated.  For  the  force  structuring  problem, 
however,  no  method  of  placing  a  bound  on  these  branches  could  be  deter¬ 
mined  and  therefore  the  method  could  not  be  used. 

Finally,  the  heuristic  algorithm  used  in  SEQOPT  was  determined 
by  the  author  and  tested  by  use  of  sample  cases. 


4.  SEQOPT  MODEL  AND  METHODS 

4.1  Method  Used  by  SEQOPT. 

SEQOPT  uses  what  may  be  described  as  a  SEQuential  Optimization 
algorithm  to  solve  the  force  structuring  problem.  It  should  be  noted  at 
the  outset  that  the  technique  used  in  this  model  has  no  connection  with 

2Ibid. 
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the  "Sequential  Unconstrained  Minimization  Technique,"  developed  by 
Fiacco  and  McCormick.  In  fact,  the  method  used  in  this  model  does  not 
in  general  produce  a  true  optimal  solution.  This  fact  will  be  illus¬ 
trated  later. 

The  sequential  optimization  technique  derives  its  name  from 
the  manner  in  which  the  force  structure  is  chosen.  Simply  put,  the 
method  chooses  the  payload  type  from  that  launch  point  which  yields  the 
maximum  target  value  damaged.  Thus,  the  method  makes  the  best  choice 
at  each  of  a  sequence  of  decision  points;  hence,  the  name  sequential 
optimization. 

Before  the  method  is  explained,  several  definitions  (in  addi¬ 
tion  to  those  in  Section  2)  are  necessary. 

V  =  ( V,  ,  V„  , . . .  :  a  vector  whose  k™1  component 

p  \  lp’  2p  TP/  th 

represents  the  value  remaining  at  the  k  target  after 
p  weapons  have  been  allocated  to  the  target  complex. 

^p  =  (^lpJ  ^2p*  *  *  •  ,c*Tp)  :  a  vector  whose  component 
represents  the  marginal  expected  damage  which  is 

th 

achievable  by  placing  one  warhead  more  upon  the  k 
target,  after  p  weapons  have  already  been  allocated. 

From  the  foregoing  definitions,  it  can  be  seen  that: 

Vo  =  (Vj,  V2,...,VT\;  and 

Dq  =  ^Vj.  Pkj,  V2  .  Pk2,...,VT  .  Pk T\. 

We  now  define  a  matrix  whose  rows  and  columns  correspond  to 
payload  type  and  launch  point,  respectively,  and  whose  elements  are  the 
marginal  expected  value  damaged  possible  from  that  launch  point  with  a 
payload  of  that  type.  Let  us  call  this  matrix  PAMVAL  (I,J),  and  now  we 
will  explain  how  its  elements  are  calculated. 


3 

Fiacco,  A.  V.,  and  McCormick,  Garth,  "The  Sequential  Unconstrained 
Minimization  Technique  for  Non-Linear  Programming,  A  Primal-Dual 
Method,"  Management  Science  10,  pp.  360-366  (1964). 
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Suppose  we  wish  to  structure  a  force  consisting  of  N  weapons , 

'til 

and  we  have  just  chosen  the  payload  type  for  the  5  weapon.  Define 

B„  =  the  number  of  warheads  which  have  been  allocated  to  the 


K 


lc  target  by  the  first  five  weapons. 


We  must  now  calculate  our  V,.  and  D_  vectors  as  follows: 

Bt  -  b 


Vk5  =  Vk 


-  t-fcj K  ]  V 


or  simplifying 


k5 


;  fork=l,...,T  and 


-  (1-PKk) 

\5={[1-(1-PKk)BK*  ]-  [-Hh]}  ■ 


or  simplifying 
^5  = 


(l-PKk)BK-(1-PKk)BK  *  ]  •  V 


;  for  k=l, . . .  ,T 


■’kS  *  PKk  (1-PKk)  K  Vk 


The  procedure  outlined  below  is  followed  in  the  calculation 
of  each  element  of  the  DAMVAL  matrix. 

(1)  Multiply  the  vector  by  A^ ,  where  i  and  j  respectively 
correspond  to  the  launch  point  and  payload  type  we  are  presently  consid¬ 
ering.  This  has  the  effect  of  yielding  no  return  for  those  targets 
which  are  not  accessible. 

(2)  Suppose  the  payload  we  are  considering  contains  m  war¬ 
heads,  then  we  do  the  following  m-times 

a.  Scan  the  vector,  selecting  the  target  whose 
marginal  return  is  the  greatest; 

b.  Since  the  value  for  the  chosen  target  has  increased 
by  1,  the  corresponding  member  in  the  D  vector  is  updated. 
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(3)  When  all  warheads  for  this  payload  type  and  launch  point 
have  been  allocated,  the  total  value  damaged  is  placed  in  the  appropriate 
position  of  the  DAMVAL  array. 


(4)  The  original  values  of  the  vector  are  replaced,  and 

the  process  is  repeated  at  (1)  if  the  DAMVAL  matrix  has  not  been  filled. 


(5)  When  the  DAMVAL  matrix  has  been  completely  filled,  the 
launch  point  and  payload  type  which  can  contribute  the  maximum  expected 
damage  is  chosen  as  a  weapon  to  be  used. 


(6)  The  Vp+j  and  D^+j  vectors  are  recalculated  and  updated, 
and  the  process  begins  again  until  the  total  number  of  weapons  have  been 


given  payloads  and  warheads  have  been  allocated  to  the  target  complex. 


An  example  here  will  help  to  clarify  the  procedure. 

4. 2  Example  Problem-Description  of  Method. 

Consider  that  in  this  problem  we  have  the  following  data: 


8  targets  (T) 

2  launch  points  (LP) 

2  missile  types  (M) 

type  1  has  2  warheads  per  missile  (dj) 
type  2  has  3  warheads  per  missile  (a0) 


Accessibility  Matrix 


Targets 

1  2  3  4  5  6  7  8 


LP1 

Missile  type  1 
Missile  type  2 

LP2 

Missile  type  1 
Missile  type  2 


10  11 
0  0  11 


0  0  0 
0  0  0 


0  10  1 
0  0  0  1 


1  1  1 
1  0  1 


Values  for  targets  1  through  8  are  as  follows: 

400.,  360.,  300.,  280.,  240.,  200.,  100, 


50.  (V-array) 


1 

0 


0 

0 
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The  probability  of  kill  is  assumed  to  be  .50  for  each  target 
(PK-array) .  We  are  also  assuming  that  any  payload  type  can  be  used  at 
either  launch  point. 

★ 


Matrix  of  Marginal  Values 

Targets 


1 

2 

3 

4 

5 

6 

7 

8 

1 

200. 

180. 

150. 

140. 

120. 

100. 

50. 

25. 

2 

100. 

90. 

75. 

70. 

60. 

50. 

25. 

12.5 

3 

50. 

45. 

37.5 

35. 

30. 

25. 

12.5 

6.25 

4 

25. 

22.5 

18.75 

17.5 

15. 

12.5 

6.25 

3.12 

Define  the  element  (i,j)  of  the  above  matrix  as: 

the  marginal  value  obtained  by  placing 
,.  ...  i  warheads  on  target  j,  over  that 

~  obtained  by  placing  i-1  warheads  on 
the  same  target 

Hence,  (i,j)  =  [(l-pk)* ”*-(l-pk)^]  .  value  of  the  j™  target. 

In  this  example  we  are  asked  to  choose  four  missiles,  so  we 
go  through  the  following  procedure  four  times: 

(1)  For  each  weapon  type  at  each  launch  point,  allocate  the 
warheads  to  that  set  of  accessible  targets  which  yield  the  greatest 
marginal  return.  The  marginal  return  is  the  sum  of  the  individual 
returns  for  this  type.  (Remember  to  multiply  the  marginal  return  array 
by  the  appropriate  accessibility  array). 

(2)  Choose  that  missile,  specified  by  launchpoint  and  type, 
which  obtained  the  maximum  marginal  return.  This  missile  is  the  choice 
for  this  loop. 

(3)  Observe  the  allocation  of  the  chosen  missile.  Suppose 
that  the  k**1  target  was  allocated  two  warheads,  then  cross  out  the  two 
top  entries  in  the  k**1  column  not  yet  removed.  The  next  value 

in  this  column  is  then  the  marginal  return  for  placing  one  more  warhead 
on  the  k^  target. 

(4)  Return  to  (1)  if  all  missiles  have  not  yet  been  chosen. 

_ 

The  model  does  not  store  this  matrix  of  marginal  values  but  computes  a 
new  entry  for  a  target  whenever  the  number  of  warheads  so  far  applied  is 
changed. 
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Choice  of  Weapons  Example 


ITERATION  1  Allocations  Marginal  Value  Obtained 


Base  1 

Type  1 

1(1,3) 

350. 

Type  2 

2(3),  1(4) 

365. 

Base  2 

Type  1 

1(2,4) 

320. 

Type  2 

2(4),  1(5) 

330. 

Hence,  choice  number  1  is  Base  1  type  2. 


ITERATION  2 

Base  1 

Type  1 

2(1) 

300.* 

Type  2 

1(3),  2(4) 

142.50 

Base  2 

Type  1 

1(2,5) 

300.* 

Type  2 

1(4),  2(5) 

250. 

It  does  not  matter  here  which  choice  is  made  because  the 
target  sets  attacked  by  the  choices  have  no  intersection.  (The  next 
choice  would  certainly  choose  that  missile  not  selected  now.)  If  any 
of  the  same  targets  were  attacked  by  both  choices,  the  first  type  would 
be  chosen,  the  previous  choices  and  alternate  would  be  written  on  tape. 
Then,  when  this  example  was  completed  using  the  first  type  as  this  choice, 
the  program  would  then  compute  following  choices  as  if  the  second  choice 
here  had  been  chosen.  In  this  manner  all  possible  choices  are  considered. 

Consider  choice  number  2  to  be  Base  1  type  1. 


ITERATION  3 


Base 

Type  1 

1(1,4) 

120. 

Type  2 

1(3),  2(4) 

142.50 

Base  2 

Type  1 

1(2,5) 

* 

300. 

Type  2 

1(4),  2(5) 

250. 

Hence,  choice  number  3  is  Base  2  type  1. 
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ITERATION  4 

Allocations 

Marginal  Value  Obtained 

Base  1 

Type  1 

1(1,4) 

120. 

Type  2 

1(3),  2(4) 

142.50 

Base  2 

Type  1 

1(2),  1(6) 

190.* 

Type  2 

1(4, 5, 7) 

180. 

Hence,  choice  number  4  is  Base  2  type  1. 


The  notation  1(1,3)  means  1  warhead  was  placed  on  targets  1  and  3. 
* 

Indicates  weapon  type  and  launch  point  selected  on  the  iteration, 
and  greatest  marginal  damage  achieved. 

Summarization  of  Outloads  and  Allocation 


Marginal 

Allocation  Expected  Damage 


Base  1 

1  Missile 

type  2 

2(3),  1(4) 

365.00 

1  Missile 

type  1 

2(1) 

300.00 

Base  2 

2  Missiles 

type  1 

1(2,5) 

300.00 

1(2,6) 

190.00 

Expected 

%  Value 

Target 

Original  Value  Damage 

Warheads  Applied 

Expected  Damage 

1 

400. 

(300.  ) 

2 

.750 

2 

360. 

(270.  ) 

2 

.750 

3 

300. 

(225.  ) 

2 

.750 

4 

280. 

(140.  ) 

1 

.500 

5 

240. 

(120.  ) 

1 

.500 

6 

200. 

(100.  ) 

1 

.500 

7 

100. 

(  0.00) 

0 

.000 

8 

50. 

(  0.00) 

0 

.000 

Total  value  of  target  base  -  1930.00 
Total  expected  damage  -  1155.00 
%  Expected  damage  to  target  base  -  59.8 

Note  here  that  the  choices  using  sequential  optimization  have  led  to  a 
feasible  solution. 
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4 . 3  Two  Phases  of  SEQOPT. 


There  are  two  distinct  phases  in  the  operation  of  this  model. 

In  Phase  I  the  constraint, 

Z  6. .  =  W.  , 

j  13  1 

is  ignored,  while  in  Phase  II  this  capacity  constraint  is  observed. 

The  ignoring  of  this  constraint  in  Phase  I  has  the  following 

effects : 

(1)  Shows  the  advantage,  if  any,  of  allowing  the  number  of 
weapons  at  any  launch  point  to  exceed  its  real  capacity; 

(2)  Shows  the  relative  worth  of  the  launch  points  by  com¬ 
paring  the  number  of  weapons  actually  placed  there  by  the  algorithm  and 

its  capacity; 

* 

(3)  Even  if  the  solutions  in  Phase  I  are  unfeasible,  they 
give  an  upper  bound  to  damages  which  could  be  generated  in  Phase  II. 

(4)  If  any  solution  determined  in  Phase  I  is  feasible,  Phase 
II  is  not  executed  since  the  results  would  duplicate  those  in  Phase  I. 

If  no  feasible  solution  is  found  in  Phase  I,  the  model  proceeds 
to  Phase  II  where  all  solutions  generated  are  feasible.  This  is  accom¬ 
plished  by  eliminating  a  launch  point  from  further  consideration  once  it 
is  filled  to  capacity. 

4 . 4  Multiple  Results. 

For  a  given  problem,  Phase  I,  II,  or  both  can  return  multiple 
results.  This  occurs  because,  at  each  decision  point  where  a  weapon  pay- 
load  and  launch  point  are  chosen,  the  choice  depends  upon  the  maximum 
value  in  the  DAMVAL  matrix.  When  this  maximum  value  is  determined,  the 
launch  point  and  payload  type  used  are  chosen.  There  is,  however,  the 
chance  that  more  than  one  element  of  this  matrix  will  contain  the  maxi¬ 
mum  value.  The  following  flow  diagram  will  make  clear  how  these  tie 
cases  are  treated: 

* 

By  unfeasible,  it  is  meant  that  some  launch  point  capacity  has  been 
exceeded;  all  other  constraints  are  handled  by  the  model. 
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The  processing  of  tie  cases  is  done  for  the  following  reasons: 

(1)  In  Phase  I,  no  tie  case  can  lead  to  an  outload  capable  of 
inflicting  damage  greater  than  other  outloads.  A  tie  could,  however, 
lead  to  a  feasible  solution,  while  it  is  possible  that  choosing  only  one 
alternative  would  not. 

(2)  In  Phase  II,  tie  cases  could  lead  to  higher  return  even 
though  all  results  would  be  feasible. 

It  should  be  noted  here  that  not  all  tie  cases  need  to  be 
examined.  No  examination  of  the  tie  is  necessary  when  the  target  sets, 
attacked  by  the  tying  weapons,  are  disjoint.  This  is  due  to  the  fact 
that  regardless  of  which  weapon  is  chosen  at  this  decision  point,  the 
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expected  damage  possible  for  the  remaining  weapon  will  not  change  at 
the  next  decision  point,  hence,  will  still  be  maximum,  and  will  be 
chosen  in  the  next  iteration. 

This  rationale  was  pointed  out  earlier  in  the  example  used  to 
describe  the  method. 

4.5  Suboptimality  of  Results. 

Since  at  present  there  is  no  feasible  method  for  finding  the 
optimal  solution  to  the  problem  while  observing  the  constraints,  we  can 
only  compare  our  results  to  an  unfeasible  solution,  which  is  an  upper 
bound,  but  not  necessarily  the  least  upper  bound  of  the  solution.  This 
upper  bound  is  computed  in  the  following  way: 

(1)  Remove  all  range  constraints  so  every  target  is  within 
range  of  every  payload  type  from  every  launch  point; 

(2)  Using  the  same  D  vector  technique  that  is  used  in  the 
algorithm,  place  one  warhead  at  a  time  on  that  target  which  yields  the 
maximum  marginal  return;  update  the  marginal  return  from  this  target 
and  continue  until  the  maximum  number  of  warheads  have  been  allocated. 

This  method  produces  the  greatest  possible  return  for  a  spe¬ 
cific  number  of  warheads,  although  in  reality  it  may  not  be  feasible. 

By  use  of  this  method,  it  is  possible  to  see  how  near  our  feasible  solu 
tion  comes  to  the  unfeasible  one  which  provides  the  upper  bound. 

At  the  outset,  it  was  felt  that  the  results  of  Phase  I,  if 

feasible,  were  an  optimum  solution;  that  is,  the  objective  function  was 

maximized.  This,  however,  was  shown  to  be  false  by  the  following 
4 

counter  example. 

Given  two  launch  points,  each  with  capacity  of  one  weapon, 
and  one  payload  type  which  has  three  warheads  per  weapon,  allocate 
these  to  a  set  of  three  targets,  with  values  V  =  (60,  80,  60),  and 
accessibility  matrix: 

4 

Weidman,  Dr.  Donald,  U.S.  Naval  Weapons  Laboratory,  Dahlgren,  Virginia 
Conference,  Naval  Weapons  Laboratory,  Subject:  SEQOPT  Model, 

November,  1970. 
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Targets 


Base _ 1 _ 2 _ 3. 

1  110 

2  011 

Also  suppose  that  the  single  shot  kill  probability  for  each 
target  is  a  constant  .5.  Then,  it  can  be  shown  that  SEQOPT  determines 
the  following  allocation: 

Target  Number  Warheads  Applied  Expected  Value  Damaged 


1 

2 

CD 

45 

(30) 

2 

3 

(3) 

70 

(70) 

3 

1 

(2) 

30 

(45) 

Total  Exp.  Value 

145 

(145) 

Damaged 

It  can  also  be  shown  that,  in  this  case,  the  optimal  alloca¬ 
tion  is  as  follows: 

Target  Number  Warheads  Applied  Expected  Value  Damaged 

12  45 

2  2  60 

3  2  45 

Total  Exp.  Value  150 

Damaged 

The  second  allocation  (see  footnote)  is  as  follows: 

Targets 
1  2  3 

Launch  Pt  1  2  1  0 

Launch  Pt  2  0  1  2  The  entries  represent 

the  allocation  of  warheads. 


The  numbers  in  parentheses  are  the  second  allocation  determined  by  the 
model.  The  return  is  exactly  the  same  due  to  the  symmetry  of  the 
problem. 
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This  counter  example  does  in  fact  show  that  the  unconstrained 
results  are  not  necessarily  optimal.  The  example  is,  however,  not  the 
specific  form  which  this  model  was  designed  to  solve.  The  example  is 
really  of  the  standard  allocation  form,  where  no  choice  of  force  struc¬ 
ture  is  possible.  It  is  doubtful  that  problems  where  a  straight  alloca¬ 
tion  is  needed  can  be  handled  adequately  by  this  model.  Dynamic  program¬ 
ing  models  and  those  based  on  heuristic  algorithms  have  been  written  for 
the  purpose  of  allocating  fixed  forces  and  should  be  utilized  for  this 
type  of  problem.1 

4.6  Other  Application  of  Sequential  Optimization  Technique. 

The  force  structuring  problem  addressed  in  this  report  is  only 
one  application  of  this  sequential  optimization  technique.  There  are 
analogous  problems  in  other  fields  such  as  economics  where  there  is  a 
limited  set  of  resources,  and  they  are  to  be  distributed  to  maximize 
return.  The  resources  present  can  be  material  or  monetary;  however, 
there  are  some  general  qualities  which  should  be  present  in  a  problem 
to  be  solved  using  this  technique: 

(1)  There  must  be  resources  of  different  types  which  could 

be  placed  at  the  points  of  distribution,  so  that  a  mixture  at  each  point 
can  be  determined; 

(2)  The  capacities  at  the  distribution  points  must  finite; 

(3)  Accessibility  of  demand  to  distribution  points  be  known; 

(4)  Associated  with  each  demand  point  there  must  be  a  return 

function,  which  follows  the  law  of  diminishing  returns.  The  variables 

present  in  this  return  function  are  amounts  of  the  resource  types  de- 

★ 

sired  at  this  demand  point; 

(5)  The  amount  of  each  resource  type  available  must  be  known. 


1Ibid. 

★ 

It  is  felt  that  the  results  of  this  method  are  nearly  optimal  due  to 
the  form  of  the  return  functions  in  that  they  follow  the  law  of 
diminishing  returns. 
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5. 


PROGRAM  SECTION 


5 . 1  Input  Variables. 

There  are  six  distinct  data  card  types  necessary  to  completely 
describe  the  variables  needed  for  a  run  of  this  model.  They  are: 

COLUMNS  FORMAT 


Card  Type  1 

1-4 

14 

5-8 

14 

9-12 

14 

13-42 

1512 

Card  Type  2 

1-12 

314 

Card  Type  3 

1-80 

8011 

Card  Type  4 

1-80 

8011 

Card  Type  5 

1-80 

10F8.0 

Card  Type  6 

1-80 

10F8.7 

VARIABLE 

NBASE  -  Number  of  launch  points 

NTYPE  -  Number  of  missile  types 
being  considered 

NTARG  -  Number  of  targets  in  target 
complex 

NRESTR(I),  1=1,  NTYPE  -  Number  of 

missiles  allowed  at  the  It*1 
launch  point 

NDESC(I) ,  1=1,  NTYPE  -  Number  of 

warheads  in  the  1^  missile 
type 

{[NALLOW  (I,J),  J=l,  NTYPE],  1=1, 

NBASE}  -  1  or  0  depending 

upon  whether  the  weapon 

type  is  allowed  at  the  1^ 
launch  point 

[NOYES  (I  ,J,K)  ,  K=1 ,  NTARG]  -  1  or  0 
depending  upon  whether  the 

K**1  target  is  within  range 

of  the  payload  type  from 

the  I*"*1  launch  point. 

VALUE (I),  1=1,  NTARG  -  Original 
target  values 

BPK(I) ,  1=1,  NTARG  -  The  single  shot 
kill  probability  associated 
with  each  target 
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5 . 2  Description  of  Important  Variables:  Not  Covered  as  Input. 


Variable  or  Array  Name 

Description 

DELDAM  (I) 

Contains  the  marginal  damage  obtainable 

til 

by  placing  one  warhead  more  on  the  I 

target  than  have  already  been  allocated 

there ; 

N COUNT  (I) 

contains  the  number  of  weapons  which 

th 

have  been  selected  for  the  I  launch 

KEEP  (I,J,K) 

point ; 

th 

the  number  of  warheads  of  the  J  weapon 

th 

type  from  the  I  launch  point  which 

th 

have  been  allocated  to  the  K  target; 

DAMAGE  (I,J) 

marginal  damage  obtainable  through  the 

th 

use  of  one  weapon  of  the  J  type  from 

th 

the  i  launch  point; 

EDSUM 

accumulator  for  the  total  expected 

damage  achievable  through  the  use  of  a 

given  out  load; 

NSUM  (I) 

the  total  number  of  warheads  allocated 

a.  -L 

to  the  I1  target; 

ED(I) 

th 

the  total  expected  damage  to  the  I 

target  for  the  allocation, 

NTARG 

EDSUM  =  l  ED  (I ) 

1=1 

5 . 3  Running  Time  and  Storage  Limitations. 

During  the  checkout  stage  of  SEQOPT,  cases  were  run  whose 
size  were  such  that  the  results  could  be  hand  verified.  These  cases 
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ranged  up  to  25  targets  and  3  launch  points,  and  3  missiles  per  launch 
point.  All  running  times  were  less  than  1  minute. 

When  SEQOPT  had  been  completely  checked  out,  a  test  case  was 
devised  which  would  compare  its  results  with  those  of  existing  tech¬ 
niques  used  to  solve  the  force  structuring  problem.  This  example  has 
200  targets,  4  launch  points,  and  3  missile  types  allowed  at  each 
launch  point.  The  running  time  for  this  example  was  10.75  minutes. 

Note  that  all  running  times  are  for  the  BRLESC  II,  a  second 
generation  machine,  and  include  compilation  time. 

It  is  felt  that,  when  converted  for  third  generation  computer 
use,  the  running  time  can  be  cut  in  half,  and  the  size  of  problems  which 
can  be  handled  will  be  greatly  increased. 

At  present  the  model  can  handle  problems  with  up  to: 

15  launch  points 
3  weapon  types 
200  targets 

The  storage  requirement  for  this  size  problem  is  only  35K,  which  leaves 
a  great  deal  of  memory  for  expansion. 
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5.4  SEQOPT  Program  Listing. 


APR ,13,71  BKLESC2  FORTRAN.  CB  15. .45 

*  SA11WC  392  107  4483  DUMBELL 

S  HA XT ( 30 ) M I  NS 

$  MAX0( 100000)LINES 

DIMENSION  NUMBON ( 15 ) , NDESC ( 3 ) , N ALLOW ( 1 5, 3 ) , NOYES ( 15, 3, 200 ) 
DIMENSION  VALUE(200),DELDAM(200),KATACK(45),IER(45),JER(45) 
DIMENSION  NSUH(200) ,  KF.EPI  15,3, 200)  ,  NKEEP(  15,  3, 200  )  ,  NCOUNT(  15) 
DIMENSION  ED( 20  0 ) , DAMAGE! 15,3) , I  I  IT ( 200 ) , TDAM < 200 > , NOTED < 15,  3 ) 
DIMENSION  FMT ( 6) , FV AR ( 3) , GMT ( 5 ) , ARRA Y ( 15, 3 ) , XRRA Y ( 200 ) 

DIMENSION  M0TED( 15,3) ,NRESTR( 15) ,BPK(200 ) 

DIMENSION  DESC(50) ,LOPNOS(50) 

DATA  (FMT( I ) , I-1,6)/10H( 1H0,7HTAP, 10HGET  , I  3, 1 0 , 4HX , 5 ( , 1 H  ,10H(I3, 
1 1H, , 1 X, 6H ) , 6X ) )/ 

DATA  (GnT(I ), 1-1,5)  /8H(1H  , 20X , 3H, 5 ( , 1H  , 8H ( I  3, 1H, , , 8H 1 X ) , 6X ) ) / 
DATA  ( F V AR ( I ), 1-1,3)  / 1H 1 , 1 H2, 1H3/ 

1000  FORMAT  (314,1512) 

1001  FORMAT  (314) 

1002  FORMAT  (8011) 

1003  FORMAT  ( 10F8.0) 

1004  FORMAT  (F3.2) 

1005  FORMAT  (1813) 

1006  FORMAT  (F10.4) 

1007  FORMAT  (16) 

1008  FORMAT  (2014) 

1009  FORMAT  ( 10(F5.4,3X) ) 

1117  FORMAT  ( 1HO, 12X, 13, 14X,F6.0, 10X, 12, 10X,F9,2, 10X,F5.3) 

102  FORMAT  ( 1 NO , 33HT0T AL  EXPECTED  DAMAGE  TO  COMPLEX-, F9 . 2 ) 

LOT  BOM-u 

LOwER-1 

NTkACK-1 

READ  10U0,  DBASE, N TYPE, NTARG, (NKESTR( I ) , I"1,NBASE) 

READ  10U1,  (NDESC( I ) , I»1,NTYPE) 

READ  1002,  (  l HALLO h ( I,J) ,J-1,NTYPE), I-l,NBASt) 

DO  6  J- 1 , NT YPE 
DO  6  I  - 1 , DBASE 

READ  1002,  ( NOYES ( 1 , J,K ) , K-l ,NT  ARG ) 

6  CONTINUE 

DO  1  I -1 ,NBASE 
1  LOTBOM-LOTBOH  +  NRESTR ( I  ) 

DO  5  I - 1 , DBASE 
5  NUmECM ( I ) -LOTBOM 

READ  1003,  ( VALUE (I),I-i,NTARG) 

READ  1009,  ( BPK ( I ) , I - 1 , NT  ARG ) 

PRINT  82 

82  FORMAT  ( ltll ,//////////, 25X, 45HBEL0W  ARE  GIVEN  THE  OPTIMAL  RETURNS 
1POSSIRLE, ,/, 

225X,40HOISREGARDING  ALL  RANGE  CONSTRAINTS,  ,FOR,/, 

325X, 45HALL  INTEGRAL  AVERAGE  NUMBERS  OF  RV  S  POSSIBLE,/, 

425X , 42HF0R  THIS  PARTICULAR  SET  OF  TARGETS.  , THESE,/, 

525X, 3 8H VALUES  ARE  UPPER  BOUNDS  FOR  NUMBERS  OF  ,/, 

625X , 42HRV  S  USED,  BUT  DOES  NOT  IMPLY  THAT  EXPECTED,/, 
?25X,49HDAMAGE  IS  MONOTON  I C ALLY  INCREASING  WITH  RV  S  USED) 

PRINT  83 

83  FORMAT  ( 1 HO, /////, 12X, 14HAVERAGE  NUMBER, 8X, 10HTOTAL  RV  S,9X, 
H5H0PTIMMM  ED  WITH,/,13X,12HRV  S/M  I SS I LE, 1 2X , 4HUSED, 1  OX , 

220HNO  RANGE  CONSTRAINTS) 

TOTKEP-O • 

DO  61  I - 1 , NT YPE 
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DE  SC ( I ) -MDESC ( I ) 

61  CONTINUE 

CALL  MAXI (DESC, NT YPE, MEMBER, BIG) 

CALL  MIn1(DESC,NTYPE,JOK,TINY) 

MAS-BIO 
MENOS-TINY 
NROKET-LOTBOM-MAS 
DO  62  1-1 f NTAR6 
PK-BPK ( I  ) 

DELPAM ( I )-PK*VALUE( I ) 

NKEEPdr  1,I)«0 

62  CONTINUE' 

KNEBO-MAS-MENOS+1 
DO  63  I « 1 , KNEBO 

LOPNOSI I ) - t  MENOS  +  I -1 )*LOTBOM 

63  CONTINUE 
LO-nLY-1 

DO  100  1 - 1 , KNEBO 
NB I GST-LOPNOS ( I ) 

DO  75  J«LOULY,NBIGST 

CAlL  MAX 1( DELOAM, NT ARG, MEM, BIG) 

NKEEP< i, l,MEM)-NKEEPt  l,l,MEMH-i 
TOTKEP-TOTKEP+BIG 
NExP-NKEF.rt  1, 1,MEM)+1 
MEXP-NEXP-1 
PK-BPK ( MEM ) 

PELDAM(mEM)«( ( ( l.-pK)**MEXP)-( ( 1 ,-PK ) **NEXP ) )*VALUE(MEM) 

75  CONTINUE 
L0wLY»N8lGST+l 

NA V-LOPNOS ( I J/LOTBOM 

PRINT  7b,  NAV,LOPNOS( I ),TOTKEP 

76  FOKMAT  ( 1HO , 1 7X, I  5, 13X, 1 5, 1 3X, F 1 0 , 2 ) 

100  CONTINUE 

4  DO  2  I-1,UTARG 
PK-BPK ( I ) 

2  DELDAM ( 1 ) -VALUE ( I )*PK 
NWKITE-I 

NREAD-2 

LIBBY-0 

NREDO-O 

DO  3  I  -1 , NBASE 
NCOUNT ( I ) -o 
DO  3  J- 1 , NT YPE 
DAMAGE  (I,J)-0, 

NOTED ( I , J ) -0 
DO  3  K-l,NTARG 
KEEP( I, J,K)»0 

3  NKEEP ( I , J, K ) -0 

10  DO  9000  KKKK-LOWER,LOTBOM 
DO  8000  KK- 1 , NBASE 

IF  (NCOUNT(KK) *LT .NUMBOM ( KK ) )  GO  TO  20 
DO  15  JJ-1,NTYPE 
DAmAGEIKK, JJ)-0, 

DO  15  I-lf NTARG 
NKEEPIKK, JJ, I ) -0 
15  CONTINUE 
GO  TO  8000 

20  DO  7000  K-l/NTYPE 

IF  (NALLOW(KK,R) .EO.l)  GO  TO  25 
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DAMAGE ( KK, K ) "0 , 

DO  21  I  ■  1 , NTARG 
21  NKtEP ( KK , K,  I )  "0 
GO  TO  7000 

25  IF  ( DAMAGE (KK,K)»Eua0.)  GO  TO  30 
GO  TO  7000 

30  DO  31  I  ■  1 , NTARG 
TEEHEE-N0YES1KK, K,  I ) 

31  TDAM( I )-DELDAM(I )*TEEHEE 
NER-NDESC(K) 

DO  6000  JJ*1,NER 

CALL  MAXI (THAN, NTARG, I TH, AMAX) 

DAMAGE  UK, K) -DAMAGE (KK,K)+TDAM ( I TH) 

NKEEPUK, K, ITH)«NKEEP(KK,K, ITH)+1 

NEwEXP-0 

DO  50  I-1,MRASE 

DO  50  J* 1 , NT YPE 

50  NEwEXP-NENEXP+KCEP< I,J, I TH ) 

NEwEXP-nF.WEXP  +  NKEEP(KK,K, ITH)+1 

MEnEXP-nEWEXP-1 

PK-BPM ITH) 

TDAMUTH)«(  (  (l,-PK)**MEWEXP)-<  ( 1 .-PK ) **NEWEXP ) ) *  V ALUE < I TH ) 

6000  CONTINUE 
7000  CONTINUE 
8000  CONTINUE 

-  - AT  THIS  POINT  HAVE  COMPUTED  FOR  ALL  TYPES  FROM  ALL  BASES  NOT  FILLED 

-  1.  DAMAGE (  BASE  TYPE) 

-  -  2.  NKEEP1BASE, TYPE, TARGETS) 

-  WILL  SELECT  BOMBER  WITH  LARGEST  DAMAGE 

CALL  MAXUM(DAMAGE,NBASE,NTYPE, AM IX, IER, JER/KONTEP) 

IF  (KONTEP.EQ. 1 )  GO  TO  8444 

•«**  ****MATERIAL  WITHIN  STARS  CONSIDERS  TIES,  WRITES  DATA  ON  TAPES 
KZT-1 

DO  8110  FJJA-2,K0NTER 

111- IER(l) 

JJl-JERl  1 ) 

1 12- IERUJJA) 

J J2« JER ( K J JA ) 

DO  8100  I - 1 , NTARG 

IF  (NKE£P( 1 1 1, JJ1, I ) ,GT*O.ANO,NKEEP( I  12, JJ2, I ) ,GT.O)  GO  TO  8101 

8100  CONTINUE 

8101  KZT-KZT+1 
IER  (  KZT  )  **  I  I  2 
JER (KZT ) ■ JJ2 

8110  CONTINUE 

IF(KZT.EQ.O)  GO  TO  8444 

-  - HENCE  KZT  WILL  BE  THE  FINAL  NUMBER  OF  CASES  PUT  ON  TAPE 

ONLY  THOSE  CASES  WHICH  ATTACKED  AT  LEAST  ONE  OF  THE  SAME  TARGETS 

-  - AS  THE  CHOSEN  BOMBER  MUST  BE  REDONE 

LIbBY«LlBRY+KZT-l 

DO  8425  KJJA-2,KZT 

DO  8400  I-1,NBASE 

IF  ( IER(KJJA) .EQ.I )  GO  TO  8340 

DO  8335  KOOKY«l,NTYPE 

8335  WRITE  ( NWR I TE, 1 005 )  ( KEEP ( I , KOOK Y, K ) , K- 1 , NT ARG ) 

GO  TO  8400 

8340  DO  8350  J«1,NTYPE 

IF  ( JER(KJJA).EQ.J)  GO  TO  8345 

WRITE  ( NWRITE, 1005)  ( KEEP (I , J, K ), K-l , NTARG ) 
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GO  TO  8350 

8345  DO  8346  K» 1 , NT ARG 

8346  1 1 IT(K)-NKEEP( I, J,K)+KEEP( I, J,K> 

WRITE  (NWRITE, 1005)  (IIIT(K),K"1,NTARG) 

8350  CONTINUE 
8400  CONTINUE 

KK-IERUJJA) 

JJ-JER(KJJA) 

DO  8415  I»1,NTARG 

IF  (NKEEP(KK, JJ, I ) .EQ.O)  GO  TO  8413 
KSUM»NKEEP(KK, JJ, I > 

DO  8410  J-l, NBASE 
DO  8410  L»l, NTYPE 
8410  KSUM«KSUM+KEEP( J,L, I) 

KSUM-KSUM+I 
MSUH-KSUM-i 
PK-BPI<(  I  ) 

TTT»( ( (l.-PK)**NSUM)-( 1  .-PK ) **KSUM) *VALUE <  I ) 

WRITE  (NWRITE, 1006)  TTT 
GO  TO  8415 

8413  WRITE  (NWRITE, 1006)  BELDAM ( I ) 

8415  CONTINUE 

DO  8423  I»1,NBASE 

IF  ( I  .EO.IER(KJJA) >  GO  TO  8421 

WRITE  (NWRITE, 1007)  NCOUNT(I) 

GO  TO  8423 

8421  NT)0«NCOUNT<  11*1 

WRITE  (NWRITE, 1007)  NTD 

8423  CONTINUE 
KTO-KKKK+1 

IF  (KTO.LE.LOTBOM)  GO  TO  8422 
KTU"999999 

8422  WRITE  ( NWR I TE, 1 007 )  KTO 
NACK»IEK (KJJA) 

MACK"JER(KJJA) 

DO  8424  I  - 1 , NBASE 
DO  8424  J«  1 , NTYPE 

IF  ( I .Eu.NACK. AND. J.EQ. MACK)  GO  TO  8420 
MOTED( I, J)«NOTED( I, J) 

GO  TO  8424 

8420  M0TED(I,J)-'I0TED(I,J)+1 

8424  CONTINUE 

WRITE  (NWRITE, 1008)  (< MOTED ( I , J) , J"l, NTYPE >, I"1 , NBASE > 

8425  CONTINUE 

****  CASE  DATA  FOR  ALL  TIES  NOW  WRITTEN  ON  TAPE  UNIT  8. 
8444  I TE" I ER (  1 ) 

JTE-JER(l) 

NOT  ED ( ITE, JTE) "NOTED ( ITE,  JTE>+1 
NCOUNT( ITE)«NCOUNT( ITE)+1 
DO  8500  I ■ 1 , NT  ARG 

IF  (NKEEP( ITE, JTE, I ) .EQ.O)  GO  TO  8500 

KEEP ( ITE, JTE, I ) -KEEP ( I TE, JTE, I ) +NKEEP ( ITE, JTE,  I  ) 

NZRA-0 

DO  8499  J-l, NBASE 
DO  8499  K-l, NTYPE 
8499  NZRA»NZKA+KEEPU,K,  I  ) 

NZRA-NZKA+1 
MZRA-NZRA-1 
PK"BPK ( I ) 
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DELDAM( I )■< ( (l.-PK)**HZRA)-(l.-PK)**NZRA)*VALUE(I ) 

8500  CONTINUE 

— — - IF  ANY  TYPES  NOT  CHOSEN  ATTACKED  ANY  TARGETS  ALSO  ATTACKED  BY 

THE  CHOSEN  TYPE,  THEIR  CONTRIBUTIONS  ARE  SET  TO  ZERO, AND  WILL 

BE  RECALCULATED 

KKuUNT-1 

DO  8600  I ■ 1 , NT ARG 

IF  (NKEEPt ITE,JTE, 1) .Ed. 0)  GO  TO  8600 
KAT  ACK ( KKOUNT ) "I 
KKOUNT-KKOUNT+1 
8600  CONTINUE 

KKOUNT-KkOUHT-1 
DO  8800  1*1, N BASE 
DO  8780  J*1 , NTYTE 
DO  8750  K«1,KK0IInT 
KTH-KATACK(K) 

IF  (NKEEPt I, J,KTH) .EQ.O)  GO  TO  8750 
DAMAGE ( I , J)-0, 

DO  e749  KB* l, NT ARG 

8749  NKEEPt I, J,KR)-0 
GO  TO  8780 

8750  CONTINUE 
8780  CONTINUE 
8800  CONTINUE 
9000  CONTINUE 

-  -  WHEN  EXIT  FROM  THIS  LOOP  HAVE  ASSIGNED  ALL  BOMBERS 

ALL  PERMANENT  ASSIGNMENTS  ARE  IN  KEEP  ARRAY 

9010  EDSUM-0 

DO  9100  I ■ 1 , NT  ARG 

NSUM ( I )-0 

DO  9050  J* 1 , NBASE 

DO  9050  K* l , NT YPE 

NSUM  (  I  )*NSU!1(  I  )  +KEEP  (  J,  K ,  I  ) 

9050  CONTINUE 

NOOO-NSU” (  I  ) 

PK*BPK(  I ) 

ED < I )■( l.-( l.-PK)**NOOO)* VALUE <  I  ) 

EDSUM*EuSIIM  +  ED<  I  ) 

9100  CONTINUE 

C— — -  NOW  BEGIN  OUTPUT  FOP  THIS  CASE 

IF  (NTRACKtEQ.2)  Go  TO  9116 
DO  9110  1*1, NBASE 

IF  (NCOyNT(I).GT.NRESTRtm  GO  TO  9115 
9110  CONTINUt 

PRINT  10110 

10110  FORMAT  < IHI,20X,45H*****THIS  SOLUTION  YIELDS  OPTIMAL  RETURN*****) 
CHECK-1. 

GO  TO  9116 

9115  PRINT  10119 

10119  FORMAT  ( 1H 1 , 20X, 40H*****TH I S  IS  AN  UNFEASIBLE  SOLUTION*****) 

9116  PRINT  10111 

10111  FORHAT  (  IH0,65HF.ELUW  ARE  THE  NUMBERS  OF  WEAPONS  OF  EACH  TYPE  CHOSE 
IN  AT  EACH  BASE) 

DO  9200  M-i, NBASE 
PRINT  10121 
10121  FORMAT  (1H0) 

PRINT  1014,  M,  NUMBOM(M) 

1014  FORHAT ( 1H  , 5HBASE  ,I2,27H  NUMBER  OF  WEAPONS  ALLOWED-, 12) 

DO  9190  1*1, NT YPE 
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PRINT  1013,  I,NOTED(M,I) 

1013  FORMAT ( 1H0, 5HTYPE  ,11, 2H  -,I2) 

9190  CONTINUE 
9200  CONTINUE 

PRINT  8062 

8062  FORMAT (1 MO/ 91HBEL0W  IS  A  LISTIN0  OF  TARGETS, NUMBER  OF  WARHEADS  APP 
1LIED,  AND  VALUE  BEFORE  AND  AFTER  ATTACK) 

PRINT  1015 

1015  FORMAT  ( 1H0, 1 5X, 6HT ARGET , 8X, 8H0R I GI NAL, 8X, 9HNUMBER  0F/8X, 
ISHEXPECTED, 8X, 13HPERCENT AGE  OF) 

PRINT  1  u  16 

1016  FORMAT  ( 30X, 5H VALUE, 6X, 16HWARHE ADS  APPLI ED, 6X , 6HDAMAGE, 9X, 

1 IShVALUE  DAMAGED) 

PRINT  10121 
DO  11111  I ■ l , NTARG 
RRINTD-EDI I )/VALUE( I ) 

PRINT  11 17, I, VALUE! I ) ,NSUM( I ) ,ED ( I ),RRINTD 
11111  CONTINUE 

PRINT  102,  EDSUM 
PRINT  1017 

1017  FORMAT  ( 1H1,73HBEL0W  IS  THE  ALLOCATION  OF  WEAPONS  TO  TARGETS,  SEPA 
1RATED  BY  BASE  AND  TYPE) 

FMT(4)-FVAR(NTYPE) 

GMT(3)-FVAR(NTYPE) 

DO  1333  1-1, NTARG 
IF  (NBASE. GT. 5)  GO  TO  1323 
NYO-NBASE 
GO  TO  1324 


1323 

1324 


1326 

1327 


1333 


«1, NTYPE), J-l.NYO) 


1400 

1410 

1420 

1425 

1430 


( (KEEP! J,K,I ), K-l, NTYPE), J-6/NY0) 
GO  TO  1333 

( (KEEP! J,K, I ),K-1,NTYPE), J-ll/NYO) 


NYU-5 

WRITE  (6,FMT)  I, ( (KEEP! J,K,I ),K 
IF  (NBASE.LE.5)  GO  TO  1333 
IF  (NBASE, LE.10)  GO  TO  1326 
NYO-10 
GO  TO  1327 
NYO-NBASE 
WRITE  I  6, GMT ) 

IF  (NBASE.LE.10) 

NYO-NBASE 
WRITE  1 6, GMT ) 

CONTINUE 
IF  (NREOO.EQ.O)  GO  TC  1500 

********«THIS  IS  READ  IN  SECTION  FROM  TAPE  NREAD 
DO  1410  I ■ 1 , NBASE 
DO  1410  J- 1 , NTYPE 

READ  (NREAD, 1005)  ( KEEP ( I , J, K ) , K-l , NT ARG ) 

DO  1420  1-1, NBASE 
DO  1420  J-l, NTYPE 
DAMAGE! I, J)-0, 

DO  1420  K-l, NTARG 

NKEEP ( I , J,K)-0 

DO  1425  1-1, NTARG 

READ  (NREAD, 1006)  DELDAM ( I ) 

DO  1430  J-l, NBASE 
READ  (NREAD, 1007) 

READ  (NREAD, 1007) 

READ  (NREAD, 1008) 

•••••••• *THIS  ENDS 

NREDO-NREDO-1 

IF  (LOWER. EQ. 999999)  GO  TO  9010 


NCOUNT! J) 

LOWER 

((N0TED(I,J),J-1,NTYPE),1' 
THE  READ  IN  SECTION 


1, NBASE) 
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GO  TO  10 

1500  IF  (LIBBY, EQ.O)  GO  TO  1600 
RE* I ND  NREAD 
REWIND  nWRITE 
NREDO-LIBRY 
LIBBY-0 
NREAD-NwRITE 

IF  (NREAD. EQ.l)  GO  TO  1502 
NWKlTE-1 
GO  TO  1505 
1502  NWRITE-2 
1505  GO  TO  1400 

1600  CONTINUE 

IF  (NTRACK.EQ.2)  GO  TO  1700 
IF  (CHECK, EO.O.)  GO  TO  1610 
PRINT  1601 

1601  FORMAT  ( IH 1 , 1 5X, 89H***** AN  OPTIMAL  FEASIBLE  SOLUTION  HAS  BEEN  FOUN 
ID,  FURTHER  CALCULATION  IS  UNNECESSARY*****) 

GO  TO  1700 
1610  LOwER  -1 

DO  1612  I»1,NRASE 
1612  NUMBOM( 1 )-NRESTR( I ) 

NTKACK-2 
REwIND  NREAD 
REWIND  NWRITE 
PRINT  1699 

1699  FORMAT  ( IH1 ,//////////, 20X, 80HN0  OPTIMAL  FEASIBLE  SOLUTION  HAS  BEE 
IN  FOUND-ALL  FOLLOWING  SOLUTIONS  ARE  FEASIBLE) 

GO  TO  4 

1700  CONTINUE 
STOP 
END 

SUBROUTINE  MAXUM( array,nrow,ncol, amix, IER, JER,KONTER) 

DIMENSION  ARRAY ( 15,3 ) , IER ( 5 ) , JER ( 5 ) 

KONTER-1 

AM1X-0. 

DO  100  l » 1 , NR0W 
DO  90  J» 1 , HCOL 

IF  (ARRAY! I,J) .LT.AHIX)  GO  TO  90 

IF  (ARRAY( I, J) .GT.AMIX)  GO  TO  80 

J  ER (KONTER ) ■ I 

JER ( KONTER ) ■ J 

KONTER-KONTER+1 

GO  TO  90 

80  AHlX» ARRAY ( I , J ) 

I ER ( 1 ) ■ 1 
JER ( 1 ) ■ J 
KONTER-2 
DO  87  K*2, 5 
IER(K)»0 
87  JER(K)-0 
90  CONTINUE 
100  CONTINUE 

KONTER-KONTER-1 

RETURN 

END 

SUBROUTINE  MAXI ( XRRAY, NSI ZE, MEMBER, THEMAX ) 

DIMENSION  XRRAY ( 200 ) 

MEMBER  ■ 1 
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THEMAX«XRRAY<1> 

DO  10  I-2,NSIZE 

IF  (XRRAY(I) .LE.THEMAX)  SO  TO  10 
THEMAX-XRRAYU  ) 

MEMBER-I 
10  continue 
RETURN 
END 

SUBROUTINE  MI Nl ( ARRAY, NS  I ZE, MEMBER, THEMI N ) 
DIMENSION  ARRAY ( 200 ) 

MEMBER-1 . 

THEHIN-ARRAY(l) 

DO  10  I«2,NSIZE 

IF  (ARRAY ( I ) .GE. THE MIN)  00  TO  10 
MEmBER-2' 

THENIN«ARRAY( I ) 

10  CONTINUE 
RETURN 
END 
LIST 

*  L I  ST ( STOP ) 
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5.5  SEQOPT  Flow  Charts . 
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5.6  Sample  Problem  Description. 


Suppose  we  have  a  target  complex  which  consists  of  ten  targets, 
with  the  following  relative  values  and  single  shot  kill  probabilities 

Targets 

1  2  3  4  ,5  6  7  8  9  10 

VALUE  300  290  260  260  210  200  190  150  100  100 

PK  .33  .42  .40  .38  .35  .30  .45  .50  .30  .37 

We  also  have  three  launch  points  each  of  which  has  a  capacity  of  four 
weapons.  We  will  consider  two  possible  weapon  types  (payloads). 

Type  I  =3  warheads/missile 
Type  II  =  2  warheads/missile 

The  accessibility  matrices  are  as  defined  below: 

Targets 


No. 

1 

Launch 

Point 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Type 

1 

0 

0 

1 

1 

0 

1 

1 

0 

0 

0 

Type 

2 

1 

1 

1 

1 

0 

1 

1 

0 

0 

1 

No. 

2 

Launch 

Point 

Type 

1 

0 

1 

0 

0 

1 

1 

1 

0 

0 

1 

Type 

2 

1 

1 

0 

0 

1 

1 

1 

1 

0 

1 

No. 

3 

Latin  ch 

Point 

Type 

1 

1 

1 

0 

0 

1 

0 

0 

1 

1 

0 

Type 

2 

1 

1 

0 

1 

1 

0 

1 

1 

1 

0 

5.7  SEQOPT  Output  For  This  Example. 

Each  iteration  of  this  model  gives  output  which  can  be  divided 
into  three  distinct  portions;  these  contain: 

(1)  The  capacity  of  each  launch  point  and  also  the  number  of 
each  weapon  type  chosen  for  each  launch  point  are  given.  (Note  that  in 
Phase  I  the  capacity  of  each  launch  point  is  equal  to  the  total  number 
of  weapons  to  be  chosen;  if  Phase  II  is  used,  the  capacities  are  reset 
to  their  true  values.) 
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(2)  The  listing  of  targets  along  with  their  original  values, 
number  of  warheads  applied  to  each  target,  expected  damage,  and  percent¬ 
age  of  target  value  damaged.  Finally,  the  total  expected  damage  for  the 
target  complex  is  given. 

(3)  The  allocation  of  warheads  to  targets,  separated  by  base 
and  type.  Turning  to  the  first  output  of  the  sample  case,  the  first 
line  is  interpreted  as  follows: 

Base  1:  types  1  and  2  had  no  warheads  delivered  to  target  1. 

Base  2:  type  1  delivered  4  warheads  to  target  1,  type  2 
delivered  0. 

Base  3:  types  1  and  2  had  no  warheads  delivered  to  target  1. 

Notice  that  Phase  I  runs  are  terminated  by  the  statement:  "NO  OPTIMAL 
FEASIBLE  SOLUTION  HAS  BEEN  FOUND  -  ALL  FOLLOWING  SOLUTIONS  ARE  FEASIBLE," 
and  all  output  following  is  from  Phase  II. 
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BELOW  ARE  GIVEN  THE  OPTIMAL  RETURNS  POSSIBLE# 
DISREGARDING  ALL  RANGE  CONSTRAINTS#  #FOR 
ALL  INTEGRAL  AVERAGE  NUMBERS  OF  RV  S  POSSIBLE 
FOR  THIS  PARTICULAR  SET  OF  TARGETS.  #THESE 
VALUES  ARE  UPPER  BOUNDS  FOR  NUMBERS  OF 
RV  S  USED,  BUT  DOES  NOT  IMPLY  THAT  EXPECTEO 
DAMAGE  IS  MONOTONICALLY  INCREASING  WITH  RV  S  USED 


AVERAGE  NUMBER 

total  RV  S 

OPTIMUM  ED  WITH 

RV  S/MISSILE 

USED 

NO  RANGE  CONSTRAINTS 

2 

24 

1441.19 

3 

36 

1704.71 

THIS  IS  An  UNFEASIBLE  SOLUTION 


bELOW  ARt  THE  IHIhBERS  OF  WEAPONS  OF  EACH  TYPE  CHOSEN  AT  EACH  BASE 

BASE  1  NUMBER  OF  WEAPONS  ALLOwEO-12 
TYPE  1  -  0 
TYPE  2  «  1 

BASE  2  Nl'IlF.EH  OF  KEATONS  ALLOWED-12 
TYPE  1  *  7 
TYPE  2  »  1 

BASE  3  n'JM&Eh  OF  WEAPONS  ALLP*ED-12 
TYPE  1  «  3 
TYPE  2-0 

BELOW  IS  A  LISTING  OF  TARGETS, NUMBER  OF  WARHEADS  APPLIED,  AND  VALUE  BEFORE  AND  AFTER  ATTACK 


TARGET 

^  r  —  T 

ORIGINAL 

value 

T  ^  —  - w 

NUMBER  OF 
warheads  applied 

EXPECTED 

DAMAGE 

PERCENTAGE  OF 
VALUE  DAMAGED 

1 

300. 

4 

239.55 

.798 

2 

290. 

4 

257*18 

,887 

3 

260. 

4 

226.30 

,870 

4 

260. 

4 

221.58 

.852 

5 

210. 

4 

172.51 

,821 

6 

200. 

4 

151,98 

,760 

7 

190. 

3 

158.39 

,834 

fl 

150. 

2 

112.50 

,750 

9 

100. 

2 

51.00 

,510 

10 

100. 

3 

75.00 

,750 

TOTAL  EXPECTED  DAMAGE  TO  COMPLEX-  1665.99 
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BELOW  IS 

TMF. 

allocation  or 

WEAPONS 

TARGET 

1 

0# 

0, 

TARGET 

2 

0, 

0, 

TARGET 

3 

0, 

Of 

TARGET 

4 

0, 

Of 

TARGET 

0, 

Of 

TARGET 

6 

0 , 

0, 

TARGET 

; 

Of 

0, 

TARGET 

3 

Of 

2/ 

TARGET 

9 

Of 

0, 

TARGET 

10 

Of 

0, 

TARGETS,  SEPARATED  BY  BASE  AND  TYPE 


4, 

0, 

Of 

Of 

4  f 

0, 

Of 

0, 

4  f 

Of 

Of 

Of 

4f 

Of 

Of 

Of 

Of 

Of 

4  f 

Of 

If 

Of 

3f 

Of 

3  f 

0, 

Of 

Of 

Of 

Of 

Of 

Of 

Cf 

2f 

Of 

Of 

If 

Of 

2f 

Of 
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* ** ** TH I S  IS  AN  UNFEASIBLE  SOLUTION***** 
bELOW  APt  THE  NUMBERS  OF  WEAPONS  OF  EACH  TYPE  CHOSEN  AT  EACH  BASE 

BASE  1  NUMBER  OF  WEAPONS  ALLOwED*l2 
TYPE  1  -  0 
TYPE  2  «  0 

bASE  2  hMMBEK  OF  WEAPONS  ALLOWED-12 
TYPE  1  «  7 
TYPE  2  «  2 

BASE  3  NUMBER  P'F  WEAPONS  Al.LOwED-12 
TYPE  1  *  3 
TYPE  2-0 


BELOW 

IS  A  LISTING  OF 

TARGETS, NUMBER  OF 

WARHEADS  APPLIED, 

AND  VALUE  BEFORE 

AND  AFTER  ATTACK 

TARGET 

ORIGINAL 

NUMBER  OF 

EXPECTED 

PERCENTAGE  OF 

VALUE 

warheads  applied 

DAMAGE 

VALUE  DAMAGED 

1 

300. 

4 

239.55 

.798 

2 

290. 

4 

257.18 

•  887 

3 

260. 

4 

226.30 

•  870 

4 

260. 

4 

221.58 

•  852 

5 

210. 

4 

172.51 

•  821 

6 

200. 

4 

151.98 

.760 

7 

190. 

3 

158.39 

•  834 

8 

150. 

2 

112.50 

.750 

9 

100. 

2 

51.00 

•  510 

10 

100. 

3 

o 

o 

• 

K. 

•  750 

TOTAL 

EXPECTED  DAMAGE 

TO  COMPLEX-  1665 

.99 
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BELOW  IS 

THE 

ALLOCATION  of 

WEAPONS  TO 

TARGET 

1 

Of 

Of 

TARGET 

2 

Of 

Of 

TARGET 

3 

Of 

Of 

TARGET 

4 

0/ 

Of 

TARGET 

5 

Of 

Of 

TARGET 

6 

Of 

Of 

TARGET 

7 

Of 

Of 

TARGET 

8 

Of 

Of 

TARGET 

9 

Of 

Of 

TARGET 

10 

Of 

Of 

TARGETS,  SEPARATED  BY  BASE  AND  TYPE 


4f 

Of 

Of 

Of 

4f 

Of 

Of 

Of 

4f 

Of 

Of 

Of 

4  f 

Of 

Of 

Of 

Of 

Of 

4f 

Of 

If 

Of 

3f 

Of 

3, 

Of 

Of 

Of 

Of 

2f 

Of 

Of 

Of 

2f 

Of 

Of 

If 

Of 

2f 

Of 
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*****this  is  AN  unfeasible  solution***** 
below  are  the  numbers  of  weapons  of  each  type  chosen  at  each  base 

BASE  t  NUMBER  Dp  WEAPONS  ALLOwED-12 
TYPE  1  -  0 
TYPE  2-0 

BASE  2  nUMHEk  OF  WEATONS  ALLOwEU«l2 
TYPE  1  «  7 
TYPE  2  «  1 

BASE  3  nUMBEK  nr  WEAPONS  allowed- 1 2 
TYPE  1  -  3 
TYPE  2  -  1 

BELOW  IS  A  LISTING  OF  TAkGETS, NUMBER  OF  WARHEADS  APPLIED,  ANU  VALUE  BEFORE  AND  AFTER  ATTACK 


TARGET 

ORIGINAL 

VALUE 

NUMBER  OF 
WARHEADS  APPLIED 

EXPECTED 

DAMAGE 

PERCENTAGE  OF 
VALUE  DAMAGED 

1 

300. 

4 

239*55 

•  796 

2 

290. 

4 

257.18 

,687 

3 

260. 

4 

226,30 

,870 

4 

260. 

4 

221,56 

.852 

5 

210. 

4 

172,51 

,621 

6 

200. 

4 

151,96 

,760 

7 

190. 

3 

156,39 

,634 

8 

150. 

2 

112,50 

,750 

9 

100, 

2 

51,00 

,510 

10 

100. 

3 

75.00 

,750 

TOTAL  EXPECTED  DAMAGE  TO  COMPLEX-  1665,99 
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bELOw  IS  THF  kLLUCATION  OF  WEAPONS  TO  TARGETS,  SEPARATED  BY  BASE  AND  TYPE 


target  1 

TARGET  2 
TARGET  3 
TARGET  4 
TARGET  5 
TARGET  6 
TARGET  7 
TARGET  8 
TARGET  9 
TARGET  10 


0,  0, 
0,  0, 
0,  0, 
0,  0, 
n,  o, 
0,  0, 
0,  0, 
0,  o, 
0,  0, 
0,  0, 


4,  0, 

4,  0, 

4,  0, 

4,  0, 

0#  0, 
1/  0, 
3,  0 , 

o,  0, 
0,  2, 
1/  0, 


0,  0, 
0,  0, 
0,  0, 
0,  0, 
4,  0, 

3,  0 1 

0,  0, 
0,  2, 
0,  0, 
2,  0, 
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****«THIS  IS  AN  UNFEASIBLE  SOLUTION***** 

BELOW  ARt  THE  NUMBERS  OF  WEAPONS  OF  EACH  TYPE  CHOSEN  AT  EACH  BASE 

base  i  number  of  weapons  allowed-12 

TYPE  1  -  0 
TYPE  2  -  1 

BASE  2  NUI1DEK  OF  WEAPONS  ALL0rED-12 
TYPE  1  ■  7 
TYPE  2-0 

BASE  3  nUMREK  OF  WEAPONS  ALLOrED-12 
TYPE  1  -  3 
TYPE  2  -  1 

BELOW  IS  A  LISTING  OF  TARGETS, NUMBtP  UF  WARHEADS  APPLIED,  AND  VALUE  BEFORE  AND  AFTER  ATTACK 


TARGET 

ORIGINAL 

VALUt 

NUMBER  OF 

warheads  applied 

EXPECTED 

DAMAGE 

PERCENTAGE  OF 
VALUE  DAMAGED 

1 

300. 

4 

239.55 

.79$ 

2 

290. 

4 

257.18 

,067 

3 

260. 

4 

226.30 

.870 

4 

260. 

4 

221.56 

.652 

5 

210. 

4 

172*51 

.621 

6 

200. 

4 

151.98 

•  760 

7 

190. 

3 

158.39 

.634 

a 

150. 

2 

112*50 

.750 

9 

100. 

2 

51*00 

,510 

10 

100. 

3 

75.00 

,750 

TOTAL  EXPECTED  DAMAGE  TO 

COMPLEX-  1665 

.99 
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BELOW  IS  THE  ALLOCATION  OF  WEAPONS  TO  TARGETS,  SEPARATE!)  BT 


BASE  AND  TYPE 


TARGET 

1 

Of 

Of 

4, 

Of 

Of 

Of 

TARGET 

2 

Of 

Of 

4, 

Of 

Of 

Of 

TARGET 

3 

Of 

Of 

4, 

Of 

Of 

Of 

TARGET 

4 

Of 

Of 

4, 

Of 

Of 

Of 

target 

5 

Of 

Of 

Of 

Of 

4f 

Of 

target 

6 

Of 

Of 

If 

Of 

3f 

Of 

TARGET 

7 

Of 

Of 

3f 

Of 

Of 

Of 

TARGET 

3 

Of 

2f 

Of 

Of 

Of 

Of 

TARGET 

9 

Of 

Of 

Of 

Of 

Of 

2f 

TARGET 

10 

Of 

Of 

If 

Of 

2f 

Of 

#****THIS  IS  AN  UNFEASIBLE  SOLUTION***** 

BELOW  ARt  THE  rillh5ERS  OF  WEAPONS  OF  EACH  TYPE  CHOSEN  AT  EACH  BASE 

BASE  1  M.JMDEK  OF  WEAPONS  ALLOaED«12 
TYPE  1  ■  0 
TYPE  2  «  0 

BASE  2  NUMBER  OF  WEAPONS  ALLOwED«l2 
TYPE  1  -  7 
TYPE  2  «  1 

BASE  3  MJMBEK  OF  WEAPONS  ALL0WED-12 
TYPE  1  «  3 
TYPE  2  «  1 

BELOW  IS  A  LISTING  OF  TARGETS, NUMBER  OF  WARHEADS  APPLIED,  AND  VALUE  BEFORE  AND  AFTER  ATTACK 


TAPGET 

original 

VALUt 

NUMBER  OF 

warheads  applied 

EXPECTED 

damage 

PERCENTAGE  OF 
VALUE  DAMAGED 

1 

300. 

4 

239.55 

,798 

2 

290. 

4 

257. IB 

,667 

3 

260. 

4 

226.30 

,670 

4 

260. 

4 

221.56 

,852 

5 

210. 

4 

172,51 

,621 

6 

200. 

4 

151,96 

,760 

7 

190. 

3 

156,39 

,634 

8 

150. 

2 

112,50 

,750 

9 

100. 

2 

51,00 

,510 

10 

100. 

3 

75.00 

,750 

TOTAL  EXPECTED  DAMAGE  TO  COMPLEX-  1663,99 
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BELOW  I 

5  THE 

ALLOCATION  OF 

WEAPONS 

TO  TARGETS, 

SEPARATED 

BY  BASE  AND  TYPE 

target 

l 

0/ 

Of 

4  f 

0, 

Of 

Of 

TARGET 

2 

Of 

Of 

4f 

0, 

Of 

Of 

TARGET 

3 

Of 

Of 

4f 

Of 

Of 

Of 

TARGET 

4 

Of 

Of 

4f 

Of 

Of 

Of 

TARGET 

5 

Of 

Of 

Cf 

Of 

4f 

Of 

TARGET 

f> 

Of 

Of 

If 

Of 

3, 

Of 

TARGET 

7 

Of 

Of 

3  f 

Of 

Of 

Of 

TARGET 

ft 

Of 

Of 

Of 

2f 

Of 

Of 

TARGET 

9 

Of 

Of 

Of 

Of 

Of 

2, 

TARGET 

10 

Of 

Of 

If 

Of 

2f 

Of 
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*****THIS  IS  AN  UNFEASIBLE  SOLUTION*#*** 

BELOW  ARt  THE  NUMBERS  OF  WEAPONS  OF  EACH  TYPE  CHOSEN  AT  EACH  BASE 

BASE  1  NUMBER  OF  WEATONS  ALLOwED-12 
TYPE  1-0 
TYPE  2-0 


BASE  2  NUMBER  OF  WEATONS  ALLOwED-12 
TYPE  1  -  7 
TYPE  2*0 


BASE  3  NUMBER  OF  wE ATOMS  ALLOrED-12 
TYPE  1  -  3 
TYPE  2-2 

bELOw  IS  A  LIFTING  OF  TARGETS, NUMBER  uF  WARHEADS  APPLIED,  AND  VALUE  BEFORE  AND  AFTER  ATTACK 

tatget 


ORIGINAL 

VALUE 


NUMBER  OF 
WARHEADS  APPLIED 


EXPECTED 

DAMAGE 


PERCENTAGE  UF 
VALUE  DAMAGED 


ft 

o 

in 


300. 
290. 
260. 
260. 
210. 
200  . 
190. 
160. 
100. 
100. 


4 

4 

4 

4 

4 

4 

3 

2 

2 

3 


239.55 

257.18 

226.30 

221.58 

172.51 

151.98 

158.39 

112.50 

51.00 

75.00 


.798 

.867 

.870 

.852 

.821 

.760 

.834 

.750 

.510 

.750 


TOTAL  EXPECTEu  DAMAGE  TO  COMPLEX- 


1855.99 
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BELOW  IS  THE  ALLOCATION  OF  WEAPONS  TO  TARGETS,  SEPARATE!!  BY  BASE  AND  TYPE 


target 

1 

0/ 

Of 

4f 

Of 

Of 

Of 

TARGET 

2 

0/ 

Of 

4f 

Of 

Of 

Of 

TARGET 

3 

0/ 

Of 

4f 

Of 

Of 

Of 

TARGET 

4 

0/ 

Of 

4f 

Of 

Of 

Of 

target 

5 

Of 

Of 

Of 

Of 

4f 

Of 

TARGET 

6 

Of 

Of 

If 

Of 

3f 

Of 

target 

7 

Of 

Of 

3f 

0, 

Of 

Of 

TARGET 

3 

Of 

Of 

Of 

Of 

Of 

2  f 

TARGET 

9 

Of 

Of 

Of 

Of 

Of 

2  f 

TARGET 

10 

Of 

Of 

If 

Of 

2  f 

Of 

NO  OPTIMAL  FEASIBLE  SOLUTION  HAS  BEEN  FOUND-ALL  FOLLOWING  SOLUTIONS  AKE  FEASI 


BELOW  ARt  THE  NUMBERS  OF  WEAPONS  OF  EACH  TYPE  CHOSEN  AT  EACH  BASE 

BASE  1  NUMBER  OF  WEAPONS  ALLOWED-  4 
TYPE  1  ■  2 
TYPE  2-2 

BASE  2  NUMBER  OF  WEAPONS  ALLOWED-  4 
TYPE  1  ■  4 
TYPE  2-0 


BASE  3  nUMREK  OF  WEAPONS  ALLOWED-  4 
TYPE  1-3 
TYPE  2  -  1 


BELOW  IS  A  LISTING  OF 

TAkGETS, NUMBER  UF 

WARHEADS  APPLIED, 

and  value  before 

AND  AFTER  ATTACK 

TARGET 

original 

NUMBER  OF 

EXPECTED 

PERCENTAGE  OF 

VALUE 

WARHEADS  APPLIED 

DAMAGE 

VALUE  DAMAGED 

1 

300. 

4 

239.55 

.798 

2 

290. 

4 

257. IS 

.887 

3 

260. 

4 

226.30 

.670 

4 

260  . 

5 

236.18 

.908 

5 

210. 

4 

172.51 

.621 

f) 

200. 

4 

151.98 

.760 

7 

190. 

3 

158.39 

.634 

A 

150. 

2 

112.50 

.750 

0 

100. 

1 

30.00 

.300 

10 

100. 

2 

60.31 

.603 

TOTAL  EXPF.CTEu  damage 

TO  COMPLEX-  1644 

.91 
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BELOW  IS  THE  ALLOCATION  OF  WEAPONS  TO  TARGETS,  SEPARATED  BY  BASE  AND  TYPE 


TARGET 

1 

0/ 

If 

3  f 

Of 

Of 

0, 

TARGET 

2 

Of 

If 

2/ 

Of 

If 

0, 

TARGET 

3 

2f 

Of 

2/ 

Of 

Of 

0, 

TARGET 

4 

3/ 

Of 

<"» 

<•  f 

Of 

Of 

0, 

TARGET 

5 

Of 

If 

0, 

Of 

3f 

0, 

TARGET 

6 

If 

Of 

If 

Of 

2f 

Of 

TARGET 

7 

0/ 

0/ 

2  f 

0/ 

If 

0, 

TARGET 

5 

0, 

If 

Of 

Of 

Of 

If 

TARGET 

9 

Of 

Of 

Cf 

Of 

Of 

If 

TARGET 

10 

0/ 

0, 

Of 

Of 

2, 

0, 

BELOW 

ARE  THE  HUNTERS  OF  WEAPONS  OF  EACH 

TYPE 

CHOSEN 

AT  EACH  BASE 

BASE 

1  nUIIHEW  of 

WE  ApO(J$ 

ALLOWED*  4 

TYPE  1 

*  3 

TYPE  2 

*  1 

BASE 

2  NUn&F.K  OF 

WEAPONS 

ALLOWED*  4 

TYPE  1 

*  4 

TYPE  2 

*  n 

BASE 

3  numrek  of 

WEAPONS 

ALLOWED-  4 

TYPE  1 

*  3 

TYPE  2 

-  1 

BELOW 

IS  A  LISTING  OF  TARGETS, NUMBER 

UF 

WARHEADS  » 

LIED, 

AND  VALUE  BEFORE 

AND  AFTER  ATTACK 

TARGET 

ORIGINAL 

NUMB  * 

OF 

expected 

PERCENTAGE  OF 

value 

WARHEADS  » 

PLIED 

DAMAGE 

VALUE  DAMAGED 

1 

300. 

4 

239.55 

.798 

p 

290. 

4 

257.18 

.887 

3 

260. 

5 

239.78 

.922 

4 

260. 

5 

236.18 

.908 

1 


5 

210. 

3 

152.33 

.725 

6 

200. 

5 

166.39 

.832 

7 

190. 

4 

172.61 

.908 

0 

150. 

2 

112.50 

.750 

9 

100. 

0 

0.00 

.000 

10 

100. 

2 

60,31 

.603 

TOTAL  EXPECTED  DAMAGE  TO  COMPLEX*  1636.83 


BELOW  IS  THE  ALLOCATION  OF 


TARGET  1  0, 
TARGET  2  0, 
TARGET  3  3, 
TARGET  4  3, 
TARGET  5  0, 
TARGET  6  2, 
TARGET  7  1, 
TARGET  8  0, 
TARGET  9  0, 
TARGET  lO  o. 


WEAPONS  TO 
If 
0, 

Of 

0, 

Of 

0, 

0, 

If 

Of 

Of 


TARGETS, 
3,  0, 

2,  0, 
2,  0, 
2,  0, 
Of  Of 
If  Of 
2,  0, 
0,  0, 
0,  0, 
Of  0, 


separated  by 
Of 
If 
Of 
Of 
3, 

2, 

1, 

Of 

Of 

2, 


BASE  AND  TYPE 

0, 

1, 

0, 

Of 

Of 

0, 

Of 

If 

Of 

Of 
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DYNAMIC  PROGRAMMING  SOLUTION  TO  FORCE  STRUCTURING  PROBLEM 


To  determine  the  dynamic  programming  solution  to  this  problem, 
it  is  first  necessary  to  formulate  the  solution  to  the  standard  alloca¬ 
tion  problem.  First,  some  additional  definitions  are  needed: 

(1)  Associated  with  each  target  there  is  a  return  function 
f^,  which  is  a  monotonically  increasing  bounded  function.  The  bound  of 
this  function  is  the  target  value  ,  hence 

Vi  =  Sup  [f.  (V.,  PK.,  N)]. 

N 

In  other  words  is  the  limiting  return  obtainable  from  the  i*"  target 
regardless  of  the  number  of  warheads  applied. 

(2)  Let  V  be  the  vector  whose  components  are  the  original 
values  of  the  targets. 

(3)  Let  A V  be  a  vector  whose  components  represent  the  remain¬ 
ing  values  of  the  targets  after  some  or  all  have  been  attacked.  Suppose 
N^  warheads  have  been  placed  on  the  i^  target,  then 


AV.  =  V. 

l  l 


f. 

l 


(V 


PK,, 


Ni) 


(4)  Let  X  be  a  two-dimensional  array  whose  elements  are 


defined  as  follows: 


X.  .  =  the  number  of  warheads  of  the  jt^1  payload 


type  placed  at  the  i^1  launch  point. 


Suppose  that  some  warheads  are  allocated.  Then  let  the  matrix 
AX  represent  the  unused  warheads. 

AX. .  =  the  number  of  warheads  of  the  payload 

^  .  th 

type  at  the  1  launch  point,  not  yet 

allocated. 


Note  that  by  specifying  the  X  matrix  the  Force  Structure  is  automatically 
determined. 
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(5)  Now  let  G  represent  the  set  of  targets  in  the  complex,  so 
G.  (t  ^ ,  *2 » *  *  *  * ^  " 

Define  a  function  H  as  follows: 


I  The  return  gained  by  distributing  the 
resources  in  X  to  any  set  of  K  targets, 
using  an  optimal  policy,  where  the 
targets  are  chosen  from  G,  having  values 
given  in  the  V  vector. 


The  recursive  relationship  is: 


Hk(X,V,G)  =  max 


max 


1,2, ...  ,t  0<XnK<Xn 
°-X12K-X12 


r  LP 

[fk<VPV  E 


M 


E  X.  .,  ) 
1=1  J=1  ljk 


+  (AX,  AV,  G  -  tR)J 


0<X 


lp,m,k!x 


LP,M 


We  know  that  Hq  (X,V,G)  =  0,  since  applying  weapons  to  no  targets  yields 
no  return . 


With  this  formulation  we  are  now  ready  to  obtain  the  solution 
to  the  force  structuring  problem.  In  this  formulation  we  have  a  two- 
dimensional  state  variable: 

State  =  (V,  Y)  , 
where  Y  is  defined  as  follows: 


y  =  (yr  y2>--->yn)  > 

’  th 

where  y^  is  the  total  number  of  weapons  of  the  i  type  which  are  avail¬ 
able  for  use.  Y  is,  therefore,  the  resource  vector. 

Our  stage  variable  will  represent  the  number  of  launch  points 
to  which  we  are  presently  assigning  weapons.  Also,  6^^ ,  will  have  the 
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same  definition  as  given  in  the  earlier  description,  that  of  defining 
the  structure  of  the  force. 


Define  our  relationship  now  as: 


Lk(Y’v) 


the  maximum  return  obtainable  from  the  set 
of  targets  for  an  optimal  choice  and  allo¬ 
cation  of  forces  from  k  launch  points 
given  a  resource  vector  Y  and  target 
values,  V°. 


The  recursion  formula  for  this  relationship  is: 


Lk(Y,V)  =  max|[HT  (yu  .  .  Nm5  V,  G)  +  Lk_1(AY  ,AV)  ]| 

0  <6kl  £  y 
0  £6k2  <_ 

0  <6km  <  y 

_  _  J  m 


It  can  be  noted  here  that  the  solution  of  the  standard  alloca¬ 
tion  problem  is  embedded  in  the  solution  of  the  force  structuring  prob¬ 
lem.  For  each  iteration  in  the  structuring  problem,  a  complete  alloca¬ 
tion  problem  must  be  solved.  For  a  problem  of  even  moderate  size  the 
computing  time  becomes  prohibitive  due  to  the  large  dimensionality  of 
the  state  variables. 

Thus,  dynamic  programming  yields  a  solution,  in  theory,  to 
the  force  structuring  problem.  However,  the  use  of  this  method  is 
impractical  due  to  time  and  memory  constraints. 
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